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We describe the results of a search for time variability of the fine structure constant a using 
absorption systems in the spectra of distant quasars. Three large optical datasets and two 2fcm/mm 
absorption systems provide four independent samples, spanning ~ 23% to 87% of the age of the 
universe. Each sample yields a smaller a in the past and the optical sample shows a 4a deviation: 
Aa/a — —0.72 ±0.18 x 10 -5 over the redshift range 0.5 < z < 3.5. We find no systematic effects 
which can explain our results. The only potentially significant systematic effects push Aa/a towards 
positive values, i.e. our results would become more significant were we to correct for them. 

PACS numbers: 98.80.Es, 06.20.Jr, 95.30.Dr, 95.30.Sf 



A common property of unified theories, applied to cos- 
mology, is that they allow space and time-dependence of 
the coupling constants ffl. Spectroscopy of gas clouds 
which intersect the sightlines to distant quasars provide 
stringent constraints on variation of the fine structure 
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constant a = #-. Observing quasars at a range of red- 
shifts provides the substantial advantage of being able to 
probe a over most of the history of the universe. 

The many-multiplet method. Variations in a would 
cause detectable shifts in the rest wavelengths of red- 
shifted UV resonance transitions seen in quasar absorp- 
tion systems. For the relativistic fine structure splitting 
in alkali-type doublets, the separation between lines is 
proportional to a 2 , so small variations in the relative sep- 
aration are proportional to a 0. The "alkali-doublet" 
(AD) method offers the advantage of being simple, but 
fails to exploit the available precision since it compares 
transitions with respect to the same ground state. In 
recent papers || Q we introduced a new technique, the 
"many-multiplet" (MM) method, which is far more sen- 
sitive than the AD method and which offers other im- 
portant advantages. The MM method allows the si- 
multaneous use of any combination of transitions from 
many multiplets, comparing transitions relative to dif- 
ferent ground-states. Simultaneously using species with 
widely differing atomic masses enhances the sensitivity 
because the difference between ground-state relativistic 
corrections can be large and even of opposite sign. The 
AD method also fails to fully exploit the available data 
since only a single doublet is analysed at a time. Using 
several different species at the same time improves the 
statistics and, importantly, provides an invaluable means 



of minimising systematic effects. 

The dependence of the observed wavenumber, lu z , on a 
is conveniently expressed as ui z — where x = 

Kfj) 2 _ V = [(fj) 4 " !]> a ° is the P re sent day value, 
and a z is the value at the absorption redshift, z. q\,q2 
are coefficients which quantify the relativistic correction 
for a particular atomic mass and electron configuration. 
These coefficients have been calculated in (|, |[ |(| using 
accurate many-body theory methods. The accuracy of 
the laboratory wavenumbers, uq, dictates the precision 
of uj z and hence the constraints on Aa/a. New high 
precision laboratory measurements of many species have 
been carried out using Fourier transform spectrographs 
specifically for the purpose of searching for varying a [fr] . 

The first application of the MM method (§ used Fell, 
Mgl and Mgll transitions in 30 absorption systems to- 
wards 17 quasars and yielded an order of magnitude 
gain over previous AD method constraints. The re- 
sults suggest a may have been smaller in the past: 
Aa/a = -1.09 ± 0.36 x 10~ 5 for 0.5 < z < 1.6, where 
Aa/a = (a z - a )/a . 

The data. In the present work, we have re-analysed 
our initial sample j|, || . Small changes in the definitions 
of the spectral fitting regions and in the selection of sys- 
tems mean we now have 28 Mg/Fell systems covering 
redshifts 0.5 < z a b s < 1.8. The Mg q coefficients are 
small compared to those for Fell, so Mg can be thought 
to act as an "anchor" against which shifts in the Fell 
lines can be measured. This large difference between the 
q coefficients enabled the dramatic sensitivity increase 
compared to the AD method. 

We include new data M, also obtained using the 
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HIRES echelle spectrograph on the Keck I telescope. The 
spectral resolution is ~ 7 km/s for the entire dataset and 
the signal-to-noise ratio per pixel is ~ 30 for most of the 
spectra. This sample is dominated by 18 damped Lyman- 
a absorption systems covering redshifts 1.8 < z a b s < 3.5 
towards 13 quasars but also includes 3 new Mg/Fell ab- 
sorption systems. Two further Keck/HIRES absorption 
systems are included (h], |ll). The redshift range is on 
average higher than the data from |4|, |) , so different tran- 
sitions are used to constrain Aa/a. The transitions used 
primarily involve multiplets of Mil, CrII and Znll. How- 
ever, other transitions (Mgl, Mgll A1II, A1III, Fell) are 
also included. Al and Si play an analogous "anchor" role 
to Mg in the lower redshift sample. 

There is an important contrast between the previous 
Mg/Fell measurements and these new ones: the Nill, 
CrII and Znll q coefficients vary not only in magnitude 
but also in sign. Some wavelengths thus shift in oppo- 
site directions for a given change in a. This, and the 
greater difference between the q coefficients (compared to 
Mg/Fell) provides a further sensitivity gain. It also di- 
lutes any possible systematic effects, especially any asso- 
ciated with wavelength calibration of the data (although 
careful tests already eliminate this as a source of signifi- 
cant error A summary of all q coefficients and all 
the wavenumbers used in our analysis, which are related 
to the same reference calibration scale, is given in tables 

i of m ra. 



A third large new optical dataset is also included in 
the present analysis. This comprises 21 SilV absorption 
doublets towards 13 quasar spectra ||. 

HI 21cm absorption lines can be compared with molec- 
ular transitions detected at mm wavelengths to constrain 
g p a (g p is the proton g-factor). We have re-analysed the 
data from [|l5| , including additional molecular absorption 
lines. This provides two new Aa/a estimates at z = 0.25 
and 0.68 (see @). 

Analysis details. The analysis methods used in the 
present work are as described in Q apart from the fol- 
lowing improvements. Aa/a is now explicitly included 
as a free parameter in a multi-parameter fit. Previ- 
ously we had varied Aa/a externally. The velocity 
width (^-parameter) of an absorption line is related to 
the FWHM of the gaseous atomic velocity distribution 
by b = FWHM/1.66, and b 2 = ^Mf + b 2 uth for an ionic 
species with mass M. The first term describes the ther- 
mal component of the line broadening at kinetic temper- 
ature, T, and the second describes a possible turbulent 
motion. T and &turb are also now included as free param- 
eters, and are not degenerate when there are > 2 species 
in a fit. Note that Aa/a and z are also not degenerate 
when there are > 2 species in a fit. We have re-analysed 
the Mgll and Fell data reported in Q using the modi- 
fied method, and the two sets of results are statistically 
indistinguishable. 

As in H, to achieve optimal precision from the data, 



all physically related parameters (z's and b's) are tied in 
the x 2 minimisation. A single z-parameter is used for 
different corresponding species. Parameter errors were 
estimated using the diagonal terms of the inverse of the 
Hessian matrix (i.e. the co-variance matrix) at the best 
fit solution. Monte Carlo simulations verified the relia- 
bility of the errors derived in this way. 

Rigorous consistency checks are imposed before a fit is 
statistically acceptable. The reduced x 2 f° r each fit must 
be ~1. Each fit is carried out in 3 different ways, first 
assuming thermal broadening (so &turb = 0), secondly 
assuming turbulent broadening (so ^jf- — 0), and thirdly 
treating feturb and T as free parameters. Variations in 
Aa/a over the 3 fits must not exceed la. Only 2 fits out 
of the optical dataset failed this test, which provides a 
simple robustness check on the derived velocity structure 
for each absorption complex. The final adopted value 
was that with the smallest reduced \ 2 (which was, as 
expected, in all but 3 cases, the third type of fit above). 

Results. We now have 72 individual estimates of 
Aa/a spanning a large redshift range, providing the most 
comprehensive constraints so far obtained. The 7 solid 
circles (annotated "many-multiplet" ) in Fig. 1. show 
the binned results for the re-analysed absorption systems 
presented in 0] and the new points based on the higher 
redshift Ni/Cr/Zn data, a total of 49 points (l3). 

The hollow triangle (annotated "alkali-doublet"), illus- 
trates the average result for the 21 SilV alkali-doublets 
]l4j . Table 1 presents a summary of the results for each 
sample. The overall deviation from Aa/a = for the 
whole optical sample is significant at the 4. la level. 

The new results for the HI 21cm/mm data are: 
Aa/a = (-0.10 ± 0.22) x 10" 5 at z = 0.25 and Aa/a = 
(-0.08 ± 0.27) x 10~ 5 at z = 0.68, assuming, without 
justification, constant g p . The error for each point in- 
cludes a component of 0.2 x 10~ 5 to allow for possible 
spatial and velocity segregation of the HI and mm ab- 
sorption. This could be due to slightly different lines of 
sight to the background quasar continuum (at such dif- 
ferent wavelengths), or differences along the same line-of 
sight, or both. A recent analysis jl7| of the same two 
absorption systems adds a systematic error ofl.7x 10~ 5 . 
Our value is derived empirically using measurements of 
the Galactic interstellar medium (see fig. 2. of Q). 

Potential systematic errors. We have carried out a 
comprehensive search for any systematic effects E^] 
which could potenatially cause the result we report. 
These include: laboratory wavelengths errors, helio- 
centric velocity variation during a quasar integration, 
isotopic saturation and abundance variation, hypcrfinc 
structure, magnetic fields, kinematic effects, wavelength 
calibration and air-vacuum wavelength conversion er- 
rors, temperature variations during the observations, line 
blending, atmospheric dispersion effects, and variations 
in the intrinsic instrumental profile. None of these are 
able to explain our result. For example, kinematic ef- 
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Sample 



Method N a 



Redshift 



Aa/a 



Fell/Mgll 

Nill/Crll/ZnII 

SilV 

21cm/mm 



MM 
MM 
AD 
radio 



28 
21 
21 
2 



0.5 < z < 1.8 
1.8 < z < 3.5 
2.0 < z < 3.0 
0.25,0.68 



-0.70 ±0.23 
-0.76 ±0.28 

-0.5 ±1.3 
-0.10 ±0.17 



TABLE I: Summary of results for 4 independent samples. 
Values of Aa/a are weighted means in units of 10 _J . MM 
and AD indicate "many-multiplet" and "alkali-doublet" . N a b s 
is the number of absorption systems in each sample. 



Redshift 



21 cm/molecular 
many— multiplet 
alkafi-doublet 




0.2 0.3 0.4 0.5 0.6 0.7 

Fractional look— back time 

FIG. 1: Aa/a vs. fractional look-back time to the Big Bang. 
The conversion between redshift and look-back time assumes 
Ho = 68 km/s/Mpc, (fi M) J2a) = (0.3,0.7), so that the age 
of the universe is 13.9 Gyr. 72 quasar absorption systems 
contribute to this binned-data plot. The hollow squares cor- 
respond to two HI 21cm and molecular absorption systems 
| |16| . Those points assume no change in g Pl so should be in- 
terpreted with caution. The 7 solid circles are binned results 
for 49 quasar absorption systems. The lower redshift points 
(below z » 1.6) are based on (Mgll/Fell) and the higher red- 
shift points on (Znll, CrII, Nill, A1III, AM, Sill) ML 28 of 
these 49 systems correspond to the sample used inM]. The 
hollow triangle represents the average over 21 quasar SilV 
absorption doublets using the alkali doublet method |14j] . 



fccts (due to velocity segregation for different species) 
could introduce a scatter in Aa/a greater than the sta- 
tistical error bars, which is not seen. Only two potentially 
significant systematic effects were identified: atmospheric 
dispersion and isotopic abundance evolution. If the spec- 
trograph slit is not parallel to the atmospheric dispersion 
direction (i.e. is not perpendicular to the horizon), dif- 
ferential dispersion will place the quasar light at different 
slit positions, depending on wavelength. In fact, this ef- 
fect turns out to push Aa/a to more positive values for 



each of the 3 optical samples. If we apply a maximum 
correction, on a case-by-case basis for the actual spectro- 
graph slit-angle, the result for the MM sample as a whole 
would become Aa/a = (-1.19 ± 0.17) x 10" 5 . 

Quasar absorption system abundances are generally 
below solar values [||, [|, so isotopic abundance ratios 
may differ from terrestrial values. Therefore, the centroid 
wavelengths for each rest-frame transition (from labora- 
tory measurements) may not be quite correct. Observa- 
tions jl^] and theoretical estimates (n]] allow us to esti- 
mate the importance of this |l^| . To do this we remove all 
weaker isotopes in all relevant species and re-fit the entire 
sample, deriving a new set of Aa/a. Again, we find that 
this effect would push Aa/a to more positive values for 
each of the 3 optical samples. If we were to apply a cor- 
rection, we would obtain Aa/a — (—0.96 ± 0.17) x 10~ 5 
for the whole MM sample. 

To summarise the above: (i) a thorough investigation 
reveals no systematic effect which can produce the our 
results, (it) applying either of the 2 significant corrections 
would enhance the significance of our results. The results 
we quote in Table 1 are not corrected for these systematic 
effects. 

Other constraints. Constraints on a variation come 
from a variety of independent sources. Laboratory mea- 
surements made over a 140 day period |^o| yield \a/a\ < 
3.7 x 10 
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yr . Another terrestrial constraint comes 

ac- 



from the Oklo natural uranium fission reactor |21 
tive ~ 1.8 x 10 9 years ago (corresponding to a "red- 
shift" of z « 0.1). Recent analyses [^2[ [2j| suggest 
Aa/a = (-0.4 ± 1.4) x 10~ 8 (although a second, sig- 
nificantly non-zero solution is also permitted adopting 
a different Sm resonance level shift). The limit above 
(favoured by p3| ) is well below our detection. The dis- 
crepancy is easily removed for a non-linear time-evolution 
in Aa/a since the quasar data probe an earlier epoch. 
Note that Fig. 1 shows that our data are consistent with 
no variation for z ,$ 1. One may also interpret the com- 
bination of the Oklo and quasar results as the absence 
of temporal variation and the existence of spatial vari- 
ation of a. Also, unlike the optical quasar data, the 
Oklo data do not constrain Aa/a directly, but constrain 
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(tq is the nuclcon-nucleon separation, 



and m-n- is the tt— meson mass). Even then, this relies 
on the unjustified assumption that the strong interaction 
and nucleon kinetic energies are constant. The Oklo re- 
sult is thus not as "clean" as the quasar results and a re- 
liable interpretation of the apparent discrepancy requires 
further work. 

Interesting limits can be obtained by comparing the 
hyperfine 21cm HI transition with optical atomic transi- 
tions in the same gas cloud. Defining X — a 2 g p m e /m p 
(m e /m p is the ratio of electron and proton masses), 
a z a b s — 1.8 gas cloud provides a limit of AX/X = 
0.7 ± 1.1 x 10~ 5 (95% confidence limit) Q. Comparison 
with our result constrains any variation of W = g p m e /m p 
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and would give a new result of AW/W = 2.1 ±0.7 x 10~ 5 
(68% limits). However, the error quoted in on AX/X 
does not include any component associated with spa- 
tial and velocity segregation, which is very likely to be 
present when comparing transitions at widely different 
frequencies, and will be important for a single measure- 
ment. The true error on AW/W is therefore probably 
significantly larger than this. 

The cosmic microwave background (CMB) probes z ~ 
1000, within ~ 10 6 years of the big bang. Future ex- 
periments p5| may reach Aa/a ~ 10~ 2 — 10~ 3 al- 
though degeneracy with any electron mass change may 
reduce this p7[ . The light element abundances constrain 
the scale lengths of additional dimensions at the time of 
primordial nucleosynthesis (z ~ 10 8 — 10 9 , a few sec- 
onds after the big bang). The 4 He yield is sensitive 
to the (uncertain) electromagnetic contribution to the 
neutron-proton mass difference p8fl . This problem is 
avoided for heavier elements and a recent analysis ]29| ] 
yields |Aa/a| < 2 x 10~ 2 . 

Interestingly, independent results are now emerging 
which support the trend in Aa/a we find. The most 
recent CMB data are consistent with a being smaller in 
the past by a few percent Also, varying speed of 

light models, plJl , are appealing because they may ex- 
plain the supernovae results for a non-zero cosmological 
constant and solve other cosmological problems (e.g. the 
horizon, flatness, monopole problems) These also 

require a smaller a in the past. We anticipate that fur- 
ther independent quasar data will provide a definitive 
check on our results. 
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